Radial sine-Gordon kinks as sources of fast breathers 



J.-G. Caputo x 'B and M.P. Soerensen 2 -0 



m 
o 



oo 



c/3 
d 



> 
in 

(N 

rn 
O 
m 



J Department of Mathematics, 
University of Arizona, 
Tucson, AZ, 85719, USA 

Department of Mathematics 
Technical University of Denmark 
DK-2800 Kgs. Lyngby, Denmark 

(Dated: March 20, 2013) 

We consider radial sine-Gordon kinks in two, three and higher dimensions. A lull two dimensional 
simulation showing that azimuthal perturbations remain small allows to reduce the problem to the 
one dimensional radial sine-Gordon equation. We solve this equation on an interval [ro,ri] and 
absorb all outgoing radiation. Before collision the kink is well described by a simple law derived 
from the conservation ol energy. In two dimensions for ro < 2, the collision disintegrates the kink 
into a last breather while for ro > 4 we obtain a kink-breather meta-stable state where breathers 
are shed at each kink "return". In three and higher dimensions d a kink-pulson state appears for 
small ro. The three states then exist as shown by a study ol the (d, ro) parameter space. On 
the application side, the kink disintegration opens the way for new types ol terahertz microwave 
generators. 

PACS numbers: 



PACS : nonlinear dynamics of solitons 05.45.Yv, mi- 
crowaves 84.40x 



I. INTRODUCTION 



The sine-Gordon equation is an important model both 
for theory and for applications. It corresponds to a classi- 
cal field with degenerated ground states (2mr) [l| . In one 
space dimension, it is integrable via the inverse-scattering 
transform and it has two main classes of solutions, the 
kink and the breather. The former is particularly inter- 
esting because it is a topological defect separating two 
regions where the solution is and 2tt. In higher di- 
mensions one can introduce the radial kink i.e. a kink 
which only depends on the radius r and this was stud- 
ied by a number of authors. Among these Christiansen 
et al @, 0], [I| have shown that such kinks with ini- 
tial velocity exhibit the return effect where they "grow" 
up to some radius and then shrink back. Note also the 
remarkable work by Geicke who described solutions of 
the radial sine-Gordon [j| @ and indicated that radial 
kinks are destroyed Q at the origin in two dimensions. 
This observation was also reported by Bogolubsky and 
Makhankhov Q . This particular phenomenon is not well 
understood. Geicke [7] reports in particular a difference 
in the collapse of the kink in two and three dimensions. 
To analyze these radial solutions one can assume that 
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the radial term is small so that the system is a perturbed 
one dimensional sine-Gordon equation. One of the an- 
alytical tools is a perturbation theory based on inverse 
scattering, see the formulation by Maslov 0. When the 
radius is small, however, this breaks down and only con- 
servation laws can be used so the analysis becomes very 
difficult. For example see the techniques used by Alfimov 
and Vazquez [l(| to study the long lived radial breathers, 
so-called "pulsons". Using this combination of analysis 
and numerics, they showed that these waves decay very 
slowly and in fact do not exist as such[i"o|. 

On the application side, the two dimensional sine- 
Gordon equation describes the electrodynamics of a 
Josephson junction between two superconducting films 
in the absence of external current and dissipation [ill ]. 
The wave part comes from Maxwell's equations and the 
sine nonlinearity from Josephson's constitutive relation. 
The variable is the phase difference (or flux) ip between 
two superconductors. In this context the kink solution, a 
"fluxon" carries a flux quantum which generates micro- 
wave radiation in the terahertz range when it collides on 
the boundary of the device. When the lateral geome- 
try of the device is reduced, the fluxon, once created, is 
"dragged" towards the narrow edge. This suggested a 
design of a particle detector [l2j and also gave rise to 
the so-called Eiffel junctions with exponential tapered 
widthfl~3|, [lj|- For these, analysis and a preliminary ex- 
perimental realization [15[ confirmed that no magnetic 
field is needed to move the kink, current alone suffices. 
The dynamics was shown to be very regular contrary to 
the standard rectangular design. Note also the analysis 
of the resonances by Jaworski (l6| . 

There is a strong link between this Eiffel design and 
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the radial sine-Gordon model as we will see below. This 
link, together with the formal studies and the applica- 
tions inspired us to undertake a numerical study of dy- 
namics of two-dimensional and higher dimensional radial 
kinks. Since the theory is difficult we relied strongly on 
numerical studies using a careful procedure. We first 
solved the two dimensional (2d) sine-Gordon equation for 
a radial kink and showed that azimuthal perturbations 
remain small. This justifies the reduction to the radial 
sine-Gordon equation. We studied this equation numeri- 
cally for a radial kink initial condition on a finite domain 
ro < r < r\, and absorbed all outgoing radiation. This 
last point is important because the radiation reflecting 
from the boundary and coming back into the compu- 
tational domain can perturb strongly the solution. We 
varied systematically ro from 10 to to see how the ra- 
dial term (d—l)(f> r /r in the Laplacian affects the collision. 
The dimension d is another parameter that we varied and 
this shows new effects. We consider the radial term as 
a perturbation of the one dimensional sine-Gordon equa- 
tion and changing ro we change the magnitude of this 
perturbation from small to very large. 

Before the collision at ro, we find that the 2d and 3d 
radial kinks are well described by a simple equation for 
the radius obtained from energy conservation. When col- 
lision occurs, the kink is always strongly affected when 
ro < 5. In two dimensions for ro < 2 it disintegrates into 
a fast breather rapidly ejected away from ro. For larger 
ro we observe a semi-stable kink-breather bound state 
which sheds fast breathers at each " return" . In all cases 
the kink decays to 0. Interestingly in three dimensions 
for ro — 5 we recover the total destruction of the kink. 
For ro < 5 all the kink energy cannot be converted into 
a single breather because the radial term is too strong 
to prevent it from escaping. Instead we observe a kink- 
pulson bound state that ejects small high-frequency (low 
energy) breathers. For d > 3 the collision can yield the 
three states as shown by a study of the (d, ro) parame- 
ter plane. This scenario explains the differences observed 
by Geike 0] and other authors for the two and three di- 
mensional kink collision. It also opens an avenue for new 
microwave devices which transform a fluxon (kink) into 
a large microwave pulse. 

The article is organized as such. In section II we illustrate 
the collapse of a sine- Gordon kink in two dimensions and 
show that there are no azimuthal effects. This justifies 
the reduction to the radial sine-Gordon equation. Sec- 
tion III is the study of its conservation laws. Using these 
we deduce a simple model for the shrinking. We examine 
in detail the radial kink collision in two, three and higher 
dimensions in section IV and characterize the emission of 
breathers. We conclude in section V and suggest a design 
for a terahertz radiation source. 



II. THE COLLAPSE OF A SINE-GORDON 
KINK IN A 2D SECTOR 

To illustrate the problem that we will consider we 
present here a 2d numerical study of the dynamics of 
a radial kink in a sector. The 2d sine-Gordon equation 
reads 



Wtt - — - Wrr - — + sin ip = 0. (1) 



As domain we consider the sector ro < r < r\ , < 8 < 8q 
shown in Fig. Q] The boundary conditions corresponding 
to no external current are homogeneous Neuman so that 
tj)0 = for 9 — 0, #o and i\) r = for r = ro, r\. 




FIG. 1: Sketch of the 2d domain for the sine-Gordon equation. 
The radial kink initial condition is shown as a dotted-dashed 
line. 



We consider the propagation of a sine-Gordon kink in- 
side such a sector. In that case the initial condition is 
given by [l[ 



0(r,t = 0) = 4atan(-^=5=), (2) 
V 1 -«o 



where Ro,uq are respectively the initial position and ve- 
locity of the kink. 
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FIG. 2: Snapshots of the evolution of a kink for the 2d sine- 
Gordon equation in a wedge. The kink is started at Ro = 12 
in a domain such that ro = 1, r\ = 20. The top panel shows 
the solution for t = 9 before collision and the bottom panel 
the solution for t = 25 after the collision. The z range for the 
top plot is [0; 2n] while for the bottom plot it is [— tt;tt]. 

We have computed the evolution of such an initial con- 
dition using the comsol finite element software [l7j . Fig. 
[5] shows two different snapshots of the evolution of a kink 
in a domain such that ro = 1, T\ — 20 assuming Rq = 10 
and uo = 0. On the top panel, showing the kink at t = 9, 
the kink is accelerated towards the narrow edge. Notice 
the absence of radiation and the characteristic overshoot. 
The bottom panel presents the solution after collision 
t = 25 and shows that the kink has disappeared and only 
a flat background <f> ks 2tt persists with some oscillations. 
Despite the violence of the collision all the energy re- 
mains in the radial n = mode and no azimuthal modes 
are excited. The space discretisation (finite elements) 
does not preserve the radial symmetry. Nevertheless az- 
imuthal perturbations do not grow. We therefore now 
look for a reduction of the model to the radial case and 
justify this approximation. 

To reduce the 2d problem it is natural to expand in 
azimuthal modes using the cosine Fourier series 

oo 

1>(r,9,t) = Y,Mr,t)Xn(9), (3) 

71=0 

where Xn{0) = cos(^f^). Plugging the expression (|3J) 
into ([TJ and projecting onto xo we obtain the evolution 

of (f> 

(b 1 f 

-0Ot t +0O rr + — = T I SUl(0o+0lXl+02X2+- • • )d0 . 

r Jo 

The integrand in the right hand side can be written as 
sin(0 o + (j>iXi + ...)= sin(0o) cos(0ixi + 02X2 + ■ • ■ ) 

+ cos(0 o ) sin(0ixi + 02X2 + •■•)■ 



The integral on the right hand side of ((4]) becomes 

^ r«» oob(^xi+ ^» +...) 

V0 Jo 

, cos(0o) [ 9o M . , , , , , . 

H 7 / of sin(0ixi + 02X2 + • ••)• 

fo Jo 

To estimate these terms we expand the cosine and sine. 
Then we find that the non zero contribution for the first 
term will yield terms of the form 

ysm(0 o ), 

and will yield cubic terms for the second integral. This 
shows that if the 4>i are small, one can assume that 

sin(0 o +0iXi+- • ■ ) ~ sin(0 o )+cos(0o)(0iXi+02X2+- • ■ ), 
so that ((4]) reduces to the radial ID sine-Gordon equation 

- (j>tt + 4>rr + — = sin(0), (5) 
r 

where the O's have been omitted for simplicity. The 
model ((5|) can be obtained for any angle 60 and in par- 
ticular for the whole two-dimensional sector. It is also 
linked to the variable width sine- Gordon equation which 
contains the term <p x w'{x)/w(x) [l3j]. The radial sine- 
Gordon corresponds to w{x) = 6ox while the Eiffel junc- 
tion is for w(x) = woe~ Xx . 



III. NUMERICAL PROCEDURE 

We now detail the numerical procedure because it is 
the basis of this work. Another reason is that the ap- 
proximate analysis based on perturbation methods is dif- 
ficult to validate a priori for small radiuses. We give it 
meaning by comparing the predictions to the numerical 
solution. We solve the radial sine-Gordon equation us- 
ing the method of lines where the space discretisation 
is done using a finite difference and the time advance 
is done using an ode solver (DOPRI5 ordinary differen- 
tial equation solver [18|]). This method is flexible and 
one can increase easily the space discretisation. Another 
time integrator we have used for comparison is the Verlet 
method [19j . The number of discretisation points for a 
typical run is 4000 and the accuracy is checked by com- 
puting the hamiltonian H (fTU|) . For all cases presented 
the relative error is smaller than 210 -5 . The boundary 
condition at ro is of the Neuman type so that there is 
perfect reflection. When ro = care must be taken be- 
cause the operator <f> r /r should be regularized because 
we have an indetermination 0/0. The way to do this is 
to invoke the limit 

0r <t>r{r,t) - 

— = 7. >r^0 4>rr(r = °> *)> ( 6 ) 
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FIG. 3: Sketch of the ID computational domain. 

so that (f) rr + r /r| r= o = 2(f>rr\r=0- 

At the instant of collision radiation is emitted from the 
kink. To avoid it re-entering the computational domain 
we introduce a "sponge layer" where waves are damped 
so that the equation becomes 



(d-1)— +sin(0) 
r 



-a(r)&, (7) 



where a(r) increases smoothly from r<j = 0.9ri to the 
edge of the domain r\ as shown in Fig. |3l This mecha- 
nism kills all radiation that travels to the right and exits 
the computational domain. The amount of energy leav- 
ing the computational domain is then computed using 
the flux relation (|T2"j) evaluated at r/ < ra- This "sponge 
layer" is better adapted for our purposes than a perfectly 
patched layer because it damps all out-going waves not 
only the ones of speed one. 

The equation (O is not integrable as in the Id case 
and there are only a finite number of conservation laws. 
These are the main analytical tool to study the solution. 
We consider them in the next section. 



IV. CONSERVATION LAWS 

The radial sine- Gordon equation in d dimensions 

<t>tt-ct>rr-(d-l)—+sm(4,)=0, (8) 
r 

in a finite domain [to,^] possesses the following energy 
conservation law. 



dH 
~dT 



(9) 



where the Hamiltonian H is 



H = f 1 r^drA + ^ + (1 - cos cj>)] = [* r^ 1 

Jro ^ 2 J ro 



drU. 
(10) 

To see this multiply © by </> t , integrate over the domain 
and integrate by parts the <j) r term. Assuming a Neuman 



boundary condition at r = ro we naturally obtain the 
flux relation for the energy 



dH 
~dt 



(11) 



By integrating this relation over time we obtain 

H(t) = H(0) + f Mn,t')Mri,t')dt' (12) 
Jo 

This enables us to compute how much energy leaves the 
computational domain at r — Tf. The energy conser- 
vation will be crucial to explain many properties of the 
solution. 

Another conservation law is related to the momentum 
II of the wave defined as 



n = / r d - Y dr$ t . 

•TO 



From the partial differential equation © we get 
dH m 



dt 



r d 1 drsai(j), 



(13) 



(14) 



which shows that even for Neuman boundary conditions 
at r — ro, r\ the momentum is not conserved. 

For a localized wave such as a kink, the integrands in 
H and II are highly localized in r. In 2d we can get 
a good approximation of the motion before collision by 
using the kink solution to the ID sine-Gordon equation 
as an ansatz 



r - R 

4>k = 4atan[exp(^= r— )], 



V(i-i? 2 )' 



(15) 



where R is the kink position. We assume that the kink 
is not too close to the boundary so the integral can be 
taken from -oo to oo. To calculate it we change the 
integration variable so that r — R = r' and write 



/+oo 
Hdr' 
-OO 



r'dr'H. 



The second integral is then because of parity. This 
gives 



H w 8 



R 



l-(R? 



(16) 



At t = we start the kink at R — Rq with uq initial 
velocity so that H = 8i?o- From that one can compute 
the velocity u — R as a function of R 



" = ±\/ l -(^) 2 (1-^) 



(17) 



Fig. 0] shows the velocity u vs the position R of a 2d 
radial kink (top panel) and the 3d radial kink (bottom 
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The first term gives 



H « 8 



I? 2 



1 - {Rf 



(20) 



The second term is a small correction of the order of the 
cube of the "width" of the kink. It is much smaller than 
the leading term (I2TJ1) . We then obtain the evolution of 
the 3d radial kink as 



(21) 



The comparison of (1211) with the numerical solution is 
also very good as shown in the bottom panel of Fig. |4] 
The differential equation for R is 



R 



^ = 

?4 U ' 



R 



(22) 



(i 



which implies a decay of R but does not have a simple 
solution. 



FIG. 4: Top panel: 2d radial kink velocity for ro = 0, 5 as a 
function of the position R. The velocity obtained from the 
numerical solution is given by the points and the analytical 
estimate (|17[) is given by the dashed line. Bottom panel: 3d 
radial kink velocity for ro = 0, 5 as a function of the position 
R. The velocity obtained from the numerical solution is given 
with the points and the analytical estimate ()21[) is given by 
the dashed line. 



panel) obtained from the numerical solution as it propa- 
gates towards ro for r = 0,1,5 and 10. The velocity is 
estimated by assuming the kink functional dependance 
(fTST). One then estimates 



Max(^) 



2u 



1-u 2 



and deduces the velocity u. On the same plot we report 
the estimate (fl~T)) . The picture shows that the agreement 
is very good even when the kink is close to ro. Assuming 
the wave keeps it's kink profile, its position R follows the 
differential equation obtained from (TTJ 



whose solution is 



R(t) — Ro cos — 2 . 



(18) 



(19) 



As expected [13J the kink travels towards the "narrow" 
end and accelerates. 

Using similar arguments we can reduce the Hamilto- 
nian for the 3d kink to 



H = R 2 Udr' 



+ 00 



r' 2 dr'n. 



V. RADIAL KINK COLLISION IN 2D AND 3D 



To understand the collision of a kink with the bound- 
ary at r = ro we have conducted extensive numerical 
studies varying systematically ro. The main result for 
both 2d, 3d and larger d is that the kink does not sur- 
vive collision in a proper way when ro is small. It decays 
after a few collisions and at each collision it emits part 
of its energy in the form of fast breathers that escape 
the radial potential. This is true for all initial conditions 
Rq > 10. This is interesting since it is a way to destroy 
a topological defect. The specifics vary from the 2d case 
to the 3d case so we will consider them separately. 



collision in 2d 



Fig. [5] shows snapshots of (f>(x, t) as a function of x for 
different times before and after the collision at ro = 0. 
The top panel shows the times t = 11, ... 36 correspond- 
ing to the kink being accelerated towards ro = 0. Notice 
the large overshoot for t = 36. The right panel of Fig. [5] 
shows the two instants t = 41 and t = 61 showing that 
very little is left of the initial kink. There is just a small 
disturbance around 2n traveling towards large r. In fact 
all the kink energy from (TlB)) H m 20 x 8 = 160 leaves 
the computational domain as shown by the energy flux 
(TT21 measured at r = 30 shown in Fig. [SJ 

The wave present in the domain after the collision is 
a fast breather. Recall that a sine-Gordon breather is 
given by 



4atan 



yl — uS 1 cos(w7(i — ur — to)) 
w cosh(-\/l — oj 2 -f(r — R — ut)) 



. (23) 
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FIG. 5: 2d radial kink collision for ro = 0. Plot of <j>(x,t) 
as a function of a; for t — 11, 16, 21, 26, 31 and 36 (top 
panel) and t = 41,61 (bottom panel). The kink is started at 
Ro = 20. The parameters are ro = 0, n = 40. 
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FIG. 6: 2d radial kink collision. Time evolution of the energy 
exiting the domain at Tf = 30, for ro = 0, 5 and 10. The kink 
is started at the same position Ro = 20 so the initial energy 
Ho ~ 160 is the same for the three cases. 



H = I67V1 - w 2 , 

so that using the same argument as for the kink we get 
for the radial case 



H b = l6R~fy/T 



(24) 



where R is the center of mass of the breather. To iden- 
tify this fast breather in the numerical solution we have 
plotted the position of its center of mass as a function of 
time on the top panel of Fig. The velocity estimated 
by the fit (dashed line) is u ~ 0.89. On the middle and 
bottom panels we have plotted the analytical expression 
of the breather (f2"3")l added to the 2tt background together 
with the numerical solution for t = 41 and t = 46. The 
parameters used for the fitting are 



0.82, t 



-3, iro = 9.2 



As can be seen the fit is very good, the error in the energy 
between the fit and the numerical solution is 10 %. Such 
a breather can "escape" the radial trap because its high 
frequency averages out the radial force. Also note that 
the radial wave equation does not support any traveling 
wave like in the ID case, so any emitted radiation has to 
be in the form of a wave packet. It is interesting to find 
breather solutions as a product of the disintegration of a 
kink. In the Id case for which the sine-Gordon equation is 
integrable, the breather and kink-antikink pairs are sep- 
arated. Here, with the radial term a connection has been 
opened between the two states. This is similar to the nu- 
merical experiments of [2l| where kink-antikink pairs are 
created out of a train of small breathers in the 4 model 
which is a perturbation of the sine-Gordon equation. 

For a larger value of the boundary ro = 5, the kink is 
reflected and looks roughly like an antikink. The snap- 
shots are shown in Fig. [5] Notice the return occurring for 
t ss 41. There is however about 20 % of energy (about 21) 
lost after the collision, as shown in the flux plot Fig. [6] 
The approximate antikink that is formed has an energy 
which is about 16 x 8 = 128 so that it "stops" at R « 16 
as shown in the bottom panel of Fig. [5J Again a breather 
is emitted, it is bound to the antikink up to t = 41 after 
which it detaches and propagates to the right. To char- 
acterize this breather we ran the simulation over a much 
larger space interval r\ — 100. The parameters of the 
breather can be calculated as follows using the energy 
loss 



where 7 is the usual Lorentz factor 

1 



7 



The energy of the breather on the infinite line is given by 



H(t = 0) - H(t = 81) w 21 = a; 16 



with xq — 30 corresponding to an instant of observation 
t = 80. This gives the following parameters 

u = 0.7, lo = 0.9995 
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FIG. 7: Characterization of the wave after 2d radial kink 
collision. The top plot shows the center of mass of the wave 
as a function of time. The middle and bottom plots show 
the numerical solution together with the breather fit (|23|l for 
t — 41 and t — 46 respectively. 

which are consistent with the energy loss (21) and the 
duration of the breather passing through the boundary 



St = width/speed = w 32. 

71* 

After the emission of this breather, the kink continues to 
oscillate and decays slowly emitting waves at each colli- 
sion with r = ro ■ For larger ro as shown in Fig. [5] for 
ro = 10 the kink decays much slower. For t — 100 its 
energy has diminished by about 10 %. At every collision 
some energy is expelled. For such value ro the radial term 
4> r /r is small and we are close to the one dimensional sit- 
uation. 




i , , , , , , 

1 20 30 40 

r 

FIG. 8: 2d radial kink collision for ro = 5. The top 
panel shows plots of (j>(x,t) as a function of x for t = 
11, 16, 21, 26, 31, 36 and 41. The bottom panel shows 
snapshots of <f>(x,t) for t = 41 and 61. The kink is started at 
Ro = 20. The parameters are ro = 5, ri = 40. 



To shed another light on the problem it is useful to 
examine the different conservation laws before and af- 
ter collision. We analyzed the difference between ro = 
and ro = 5, by computing the momentum II, the right 
hand side S in the flux of the momentum (1141) and the 
front position R. The latter is defined as the maximum 
of 4> x . The time evolution of the three quantities II, S 
and R are shown in Fig. [5] The top panel corresponds 
to r = and the bottom panel to r = 5. For r = 
the momentum II goes through for R = (collision 
instant). There S > so that II will keep decreasing 
and oscillate around zero, indicating that the soliton is 
destroyed. For ro = 5 shown in the bottom panel, II > 
at collision and S > so that II remains positive for 
t > 20. After that instant S starts to oscillate with a 
period about 6. S is largely positive so that II < on 
average. We then reach II = so that the kink stops 
at t « 41 for which R = 15. This is the return effect. 
The total energy is conserved but gets distributed differ- 
ently between the different components, the kinetic term 
e< = f r * O.50jX dx, the gradient term e x = 0.5x(j) x dx 
and the potential term e p = /^(l — cos<p)x dx. Initially 
the kink has velocity so that e* = 0, all the energy is 
concentrated in e x and e p . For ro = on the top panel of 
Fig. 1101 the kinetic term et increases from to its max- 
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FIG. 9: 2d radial kink collision. Time evolution of the 
rescaled momentum II/2, the kink position R, and the source 
term S of the momentum equation for ro = (top panel) and 
ro = 5 (bottom panel). 



imum at the collision and then remains about constant. 
The potential term e p decreases from its maximum value 
at t = and stabilizes around half its value. The behav- 
ior is different for ro = 10 shown on the bottom panel of 
Fig. [TD1 for which the collision is almost elastic. There 
at the instant of collision, the kinetic energy reaches its 
maximum, the total energy and the potential energy and 
gradient energies are almost zero. After collision both 
recover their initial values. 
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FIG. 10: 2d radial kink collision: Time evolution of the energy 
components et = f^ 1 0.5x<t>tdx, e x = 0.5xcj>^.dx and e p = 
(1 — cos<j>)xdx for a kink. On the top panel ro = while 
on the bottom panel panel ro = 10. The kink is started at 
the same position Ro — 20. 





















^cos , 







B. Collision in 3d and higher d 

We now consider the collision of a kink in 3d and higher 
d to see if there are particular situations. The general re- 
sult that radiation is emitted as fast breathers out of the 
computational domain still holds. As an example con- 
sider the flux of the energy shown for the 3d case in Fig. 
IITI for ro = 0, 5 and 10. Interestingly the case ro = 5 
is similar to the 2d case for ro = 0. The kink is entirely 
transformed into a fast breather that exits the computa- 
tional domain. This is shown in the series of snapshots 
in Fig. [T^] The fast breather is clearly seen traveling to 
the right at time t = 31. For r — there is a bound 
state kink/breather so that the kink "sheds" a breather 
at every collision with the boundary and decays. Fig. 
[TBI shows the successive snapshots of the solution in this 
case. The solution for t = 31 is clearly a combination of a 
kink and a breather. For such a large value of d and such 
a small ro the radial term (d — l)<fi r /r is very strong and 
prevents a low frequency breather from escaping. Only 
breathers of frequency u> ks 1 can escape and these have 
fairly low energy. Such a breather will be "shed" from 
the kink as it reaches its return point, around R = 5. 
The energy and momentum behave in a very similar way 




50 100 

t 

FIG. 11: Time evolution of the energy flux for three different 
values ro = 0, 5 and 10. The kink is started at the same 
position Ro — 20 so the initial energy Ho ~ 3200 is the same 
for the three cases. 



as in the 2d case so we do not present them. 

To confirm these findings we conducted two simula- 
tions with d = 5, for ro = and 5 and r\ =40. The flux 
of energy exiting the domain is shown in Fig. 1141 It shows 
two breathers being emitted respectively at t = 55 and 
t = 100 for ro = 0. The (d, ro) parameter plane is shown 
in Fig. [15] It shows the coexistence of the three states, 
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FIG. 12: 3d radial kink collision for ro = 5: plot of 4>{x,t) 
as a function of x for t = 11, 16, 21, 23, 26 and 31. The 
snapshots after the collision are indicated in dashed line (blue 
online). The kink is started at Ro = 20. The parameters are 
ro = 5,n = 40. 



FIG. 15: Parameter plane (d, ro) showing the different possi- 
ble outcomes for a kink collision, a fast breather (B), a kink- 
pulson metastable state (KP) and a kink-breather metastable 
state (KB). 
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FIG. 13: 3d radial kink collision for ro = 0: plot of <j>(x,t) 
as a function of x for t — 11, 16, 21, 23, 26 and 31. The 
field for the snapshot t — 23 has been divided by 5 to fit in 
the plot. The kink is started at Ro = 20. The parameters are 
r — 0, n = 40. 
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FIG. 14: Energy exiting the computational domain [ro : 100] 
for d = 5 and ro = in continuous line (red online) and ro = 5 
in dashed line (blue online). The total energy is 1.295743 
xlO 6 . 



the fast breather (B), the kink-breather metastable (KB) 
and the kink-pulson metastable state (KP) Interestingly 
the last one cannot be seen for d — 2. 



VI. CONCLUSION 



Motivated by theory and applications, we studied ra- 
dial sine- Gordon kinks in two and higher dimensions. A 
full two dimensional simulation showed that azimuthal 
perturbations remain small. We therefore reduced the 
problem to the one dimensional radial sine-Gordon equa- 
tion which we solve on an interval [ro, r\\. Before collision 
the kink is well described by a simple law derived from 
the conservation of energy. In two dimensions the colli- 
sion of the kink with the boundary ro will result in a fast 
breather for small ro and in a kink-breather metastable 
state for larger ro. In the latter, the kink sheds at each 
"return" a large part of its energy into bursts which are 
breather solutions. We have characterized these waves in 
terms of their energy, frequency and velocity. In three 
and higher dimensions and small ro we observe a kink 
pulson bound state. The three states exist in the (d, ro) 
parameter space. This study shows that radial perturba- 
tion opens a channel between the kink solutions and the 
breather solutions. This is particularly interesting be- 
cause in one dimension these are completely separated. 
This additional term provides therefore a mechanism to 
"destroy" these topological defects and extract the en- 
ergy they contain. 

In view of applications to 2d Josephson junctions this 
could be very useful to generate Terahertz radiation. At 
this time output from these devices is low. Here for small 
ro all the kink energy is converted into radiation. This 
suggests the design of a new device based on window 
Josephson junctions. A sketch of the device is shown in 
Fig. 1161 The top panel shows a top view of the junction 
together with the radio- frequency detector D. Notice the 
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FIG. 16: Sketch of a microwave generator based on a radial 
Josephson window junction. 

current input similar as in [l5j with its passive region sep- 
arating the electrodes. The bottom panel shows a side 
view of the system with the oxide layer O separating the 
two superconducting films. As the current I is increased 
a train of fluxons is formed. These reflect into the narrow 



end r = ro and fast breathers are formed which consist 
in bursts of microwave. Since all the kink energy is con- 
verted into microwaves, we expect this system to generate 
much more radiation than a standard flux-flow. 
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